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The ^(e, e'7r + )n cross section was measured at four-momentum transfers of Q 2 = 1.60 and 2.45 
GeV 2 at an invariant mass of the photon nucleon system of W=2.22 GeV. The charged pion form 
factor (F n ) was extracted from the data by comparing the separated longitudinal pion electropro- 
duction cross section to a Regge model prediction in which F n is a free parameter. The results 
indicate that the pion form factor deviates from the charge-radius constrained monopole form at 
these values of Q 2 by one sigma, but is still far from its perturbative Quantum Chromo-Dynamics 
prediction. 



A fundamental challenge in nuclear physics is the de- 
scription of hadrons in terms of the constituents of the 
underlying theory of strong interactions, quarks and glu- 
ons. Properties such as the total charge and magnetic 
moments are well described in a constituent quark frame- 
work, which effectively takes into account quark-gluon 
interactions. However, charge and current distributions, 
which are more sensitive to the underlying dynamic pro- 
cesses, are not well described. 

Hadronic form factors provide important information 
about hadronic structure. The coupling of a virtual pho- 
ton to structureless particles is completely determined by 
their charge and magnetic moments. However, for com- 
posite particles one must account for the internal struc- 
ture, which is accomplished by momentum transfer de- 
pendent functions. Examples of these functions are the 
electromagnetic form factors, which describe the distri- 
bution of charge and current. 

One of the simplest hadronic systems available for 
study is the pion, whose valence structure is a bound 



state of a quark and an antiquark. The electromagnetic 
structure of a spinless particle such as the pion is param- 
eterized by a single form factor. The pion charge form 
factor, F„, can be calculated in perturbative Quantum 
Chromo-Dynamics (pQCD) in the limit of very large val- 
ues of four-momentum transfer squared, Q 2 

F v (Q 2 ) = Sir ^ (Q 2 - oo), (1) 

where a s is the strong coupling constant. The normaliza- 
tion is hxed by the pion decay constant, f n , which is de- 
termined from the weak decay of the pion (tt — > fx + v^). 
At low values of Q 2 , Vector Meson Dominance models 
provide a reasonable description of F„. Due to the pion's 
simple qq valence structure, the transition from "soft" 
(nonperturbative) to "hard" (perturbative) physics is ex- 
pected to occur at significantly lower values of Q 2 for F^ 
than for the nucleon form factors [2| ■ 

The form factor of the pion is well determined up to 
Q 2 =0.28 GeV 2 by elastic tt — e scattering experiments 
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|3j, from which the charge radius has been extracted. 
Extending the measurement of to larger values of Q 2 
requires the use of pion electroproduction from a nucleon 
target. The pion exchange (i-pole) process, in which a 
virtual photon couples to a virtual pion inside the nu- 
cleon, dominates the longitudinal pion electroproduction 
cross section, ox, at small values of the Mandelstam vari- 
able t. There ol exhibits a characteristic ^-dependence 
and is proportional to F 2 . 

Experimental values of F n havepreviously been de- 
termined at CEA and Cornell 0, [| , DESY [10, and 
recently at Jefferson Lab [j| . Most of the high Q 2 data 
have come from experiments at Cornell covering a range 
of values in Q 2 between 0.28 and 9.77 GeV 2 . In these ex- 
periments Fn was extracted from the longitudinal cross 
sections, which were isolated by subtracting a model of 
the transverse contribution from the unseparated cross 
sections. Pion electroproduction data were also obtained 
at DESY for a value of Q 2 of 0.7 GeV 2 , at an invariant 
mass of the photon nucleon system of W=2.19 GeV, and 
longitudinal and transverse cross sections were extracted 
using the Rosenbluth separation method. In 1997, Jeffer- 
son Lab provided the first high precision pion electropro- 
duction data for F^ for values of Q 2 between Q 2 =0.6 and 
1.6 GeV 2 at W=1.95 GeV II. For an updated analysis 
of these data see reference |9j. These data give a precise 
determination of 07, with a significant improvement in 
the experimental uncertainty. The results presented here 
extend the Q 2 range to 2.45 GeV 2 and address questions 
of model dependence in the extraction of F n . 

The experiment described here was carried out in Hall 
C at Thomas Jefferson National Accelerator Facility (Jef- 
ferson Lab). Pion electroproduction cross sections were 
measured from hydrogen and deuterium targets. The 
data were taken at two beam energies for each of the 
two values of Q 2 at W=2.22 GeV. Charged pions were 
detected in the High Momentum Spectrometer (HMS), 
while the scattered electrons were detected in the Short 
Orbit Spectrometer (SOS). Both spectrometers include 
two drift chambers for track reconstruction and scintil- 
lator arrays for triggering. A detailed description of the 
Jefferson Lab Hall C spectrometers can be found in ref- 
erence fioj. 

In order to select electrons in the SOS, a gas Cerenkov 
detector containing Freon-12 at atmospheric pressure was 
used in combination with a lead-glass calorimeter. Posi- 
tively charged pions were identified in the HMS using an 
aerogel Cerenkov detector with refractive index of 1.03 
[nj . In the case of pion production at negative polar- 
ity, electrons were rejected using a gas Cerenkov detector 
containing C4F10 at 0.47 atm. Any remaining contami- 
nation from real electron-proton coincidences was elimi- 
nated with a coincidence time cut of ±1 ns. Background 
from alumininum target cell walls (2-4% of the yield) 
and random coincidences (~1%) were subtracted from 
the charge normalized yields. The exclusive neutron fi- 
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FIG. 1: Representative example of the measured cross sec- 
tions, as a function of cj> at Q 2 =1.6 GeV 2 for two values 
of e. The curves shown represent the model cross section used 
in the Monte Carlo simulation. 



nal state was selected with a cut on the reconstructed 
missing mass. The relevant electroproduction kinematic 
variables Q 2 ,W and t were reconstructed from the mea- 
sured spectrometer quantities. Experimental yields were 
calculated after correcting for several inefficiencies, the 
dominant sources being particle tracking efficiency (3- 
4%), pion absorption (4.8%), and computer dead time 
(1-11%). The net uncertainty in these corrections is dom- 
inated by the uncertainty in the absorption of the pions 
(-2%). 

The unpolarized pion electroproduction cross section 
can be written as the product of a virtual photon flux 
factor and a virtual photon cross section, 



d 5 <j 



dQ e dE'dfl n 



j {t,4> —> 0^) 



d 2 a 
dtd<j> ' 



(2) 



where J (t, <fi — > f2 T ) is the Jacobian of the transfor- 
mation from dtdcj) to dCl^, <f> is the azimuthal angle 
between the scattering and the reaction plane, and 

^v =1 ^TT~Q5TZi W 2M 1 1S tne virtual photon flux fac- 
tor. The virtual photon cross section can be expressed 
in terms of contributions from transversely and longitu- 
dinally polarized photons, 



2tt- 



dtdcj) 



da L 
~~dt 



Here, 



1 



d<TT 

~dT 

do~TT n , 

e — ; — cosZcp. 
dt 

2^1 



—-da LT 

e) : COSfi 

; dt 



(3) 



tan 



2 



is the virtual photon po- 
larization, where q 2 is the square of the three- momentum 
transferred to the nucleon and 9 is the electron scatter- 
ing angle. The individual components in equation|21were 
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determined from a simultaneous fit to the <fi dependence 

.2 

of the measured cross sections, 4r?j, at two values of e. 
A representative example as a function of <f> is shown in 
Figure Q 

The separated cross sections are determined at fixed 
values of W, Q 2 and t, common for both high and low 
values of e. However, the acceptance covers a range in 
these quantities, thus the measured yields represent an 
average over that range. Note that each i-bin has a differ- 
ent average value of Q 2 , W. In order to minimize errors 
resulting from averaging, the experimental cross sections 
were calculated by comparing the experimental yields to 
a Monte Carlo simulation of the experiment. To account 
for variations of the cross section across the acceptance 
the simulation uses a 1 H(e,e / 7r + )n model based on pion 
electroproduction data. In addition, the Monte Carlo in- 
cludes a detailed description of the spectrometers, mul- 
tiple scattering, ionization energy loss, pion decay, and 
radiative processes. The separated cross sections, <tl and 
<Tt, are shown in Figure 

The uncertainty in the separated cross sections has 
both statistical and systematic sources. The statistical 
uncertainty in or + €<tl ranges between 1 and 2%. Sys- 
tematic uncertainties that are uncorrelated between high 
and low e points are amplified by a factor of 1/Ae in 
the L-T separation. Correlated systematic uncertainties 
propagate directly into the separated cross sections. Un- 
certainties in the scattering kinematics and beam energy 
were parameterized using data from the over-constrained 
1 H(e, e'p) reaction. Beam energy and spectrometer mo- 
menta were determined to 0.1% while the spectrometer 
angles were determined to ~0.5 mrad. The spectrometer 
acceptance was verified to better than 2% by compar- 
ing e — p elastic scattering data to a global parameter- 
ization jl2j. The uncorrelated systematic uncertainty is 
dominated by acceptance (0.6-1.1%) resulting in a total 
uncorrelated uncertainty of 0.9 to 1.2%. The correlated 
systematic uncertainty is mainly due to radiative correc- 
tions (2%), pion absorption (2%), and pion decay (1%) 
resulting in a total correlated uncertainty of 3.5%. A 
third category of systematic uncertainties consists of un- 
certainties that differ in size between e points, but may 
influence the ^-dependence at a fixed value of e in a corre- 
lated way. The "t-correlated" uncertainty is dominated 
by acceptance (0.6%), kinematics (0.8-1.1%) and model 
dependence (1.1-1.3%) resulting in a partially correlated 
uncertainty of 1.8 and 1.9%. 

In order to determine F % , the experimental results for 
Oh are compared to a Regge model calculation by Van- 
derhaeghen, Guidal and Laget (VGL) 0. In this ap- 
proach, pion electroproduction is described as the ex- 
change of Regge trajectories for 7r- and p-like particles. 
Since most model parameters are fixed by pion photopro- 
duction data, F n and the np^ transition form factor are 
the only free parameters. Both form factors are parame- 
terized by a monopole form, [1 + Q 2 / A 2 } -1 , but the cutoff 
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FIG. 2: Separated cross sections, ol and ut at central values 
of Q 2 =1.60 (2.45) GeV 2 . Note that the average values of W 
and Q 2 are different for each —t-bin. The error bar indicates 
statistical and uncorrelated systematic uncertainty in both e 
and —t combined in quadrature. The error band denote the 
correlated part of the systematic uncertainty by which all data 
points move collectively. The curves denote VGL Regge calcu- 
lations for a l (solid line) and ax (dashed line) with values of 
Al=0.513 (0.491) GeV 2 and A 2 p =l.l GeV 2 . The discontinu- 
ities in the at curve result from the different average values 
of W and Q 2 for the various t-bins. 



Q 2 (GeV^ 




1.60 
2.45 


0.243±0.012 
0.167±0.010 



TABLE I: Extracted values for F n at a value of W =2.22 
GeV. The error on F w combines statistical and experimental 
systematic uncertainties in quadrature. 



parameter, A 2 , is not as well constrained as the pion cut- 
off parameter, A 2 . Varying A 2 between 0.6 and 2.1 GeV 2 
changes cr T by 13% (30%) at Q 2 of 1.60 (2.45), but has 
little influence on err,. Thus, F v can be determined in a 
one parameter fit from a comparison of the longitudinal 
experimental cross section to the one predicted from the 
Regge model. 

A comparison of our data to the VGL prediction is 
shown in Figure El The i-dependence of the longitudinal 
cross section is well described at both central values of 
Q 2 . However, the transverse cross section is underpre- 
dicted systematically. The value of F n was determined 
from a least squares fit of the Regge model prediction to 
the data, and the resulting values are shown in Table [I] 

The extraction of F^ from ol relies on the domi- 
nance of the pion exchange term. To test the pole dom- 
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FIG. 3: Pion form factor as extracted in this work. Also 
shown are e — n elastic data from CERN, and earlier pion elec- 
troproduction data from DESY and Jefferson Lab. The ear- 
lier Jefferson Lab data are taken from reference 0/. The data 
point at Q 2 = 1.60 GeV 2 from [JJ/ has been shifted from its 
central value for visual representation. The curves are from a 
Dyson-Schwmger equation (solid, fig). QCD sum rules (dot- 
ted, ]lAf). dispersion relations with QCD constraint ( dashed, 
and from a pQCD calculation (dashed- dotted, \l&f). 



inance the longitudinal ir~ /ir + ratios in 2 H were exam- 
ined. Since the pole term is purely isovector this ratio is 
expected to be close to unity and a significant deviation 
from unity would indicate the presence of an isoscalar 
background. The preliminary analysis of the longitudi- 
nal ir~ /tt + ratios is consistent with unity. 

In Figure |3 our results are shown along with re- 
sults from CERN, DESY, earlier Jefferson Lab data, and 
some representative calculations. Comparing the result 
at Q 2 — 1.60 GeV 2 to the earlier Jefferson Lab data 
point at a lower value of W allows for a direct test of the 
theoretical model dependence. A higher value of W al- 
lows for a measurement at smaller values of — t, at closer 
proximity to the pion pole. The data are consistent with 
the previous Jefferson Lab F n measurement at a value of 
Q 2 — 1.60 GeV 2 and suggest a small model uncertainty 
due to fitting the VGL model to the data. The data in- 
dicate a one sigma deviation from a monopolc form fac- 
tor that yields the measured charge radius. That form 
factor is up to Q 2 =2.5 GeV 2 indistinguishable from the 
solid curve in Figure |21 Various models provide a good 
description of the measured values for F T up to Q 2 =1.60 
GeV 2 . The data are well described by the calculation of 
Nesterenko and Radyushkin (l4|. in which a QCD sum 
rule framework for the soft contribution to F n as well as 
an asymptotically dominant hard gluon exchange term 
is used. The dispersion relation calculation by Geshken- 



bein [15j also agrees well with the data. The data are 
also reasonably well described by the Dyson-Schwinger 
calculation by Maris and Tandy, which is based on the 
Bethe-Salpeter equation with dressed quark and gluon 
propagators. All parameters in the latter calculation are 
determined without the use of F^ data 0,0]- Perturba- 
tive QCD calculations of which one is shown in Figure |3| 
give values of Q 2 F 7T around 0.10 GeV 2 in the region of 
our measurements. 

In summary, we have measured separated 1 H(e,e / 7r + )n 
cross sections at values of Q 2 =1.60 and 2.45 GeV 2 at 
IV=2.22 GeV. The charged pion form factor was ex- 
tracted from the separated longitudinal cross section us- 
ing a Regge model. The data are consistent with the 
previous Jefferson Lab result at Q 2 = 1.60 GeV 2 . The 
data deviate by one sigma from a monopole form factor 
obeying the measured charge radius, but are still far from 
the values expected from pQCD calculations. 

This work was supported in part by the U.S. Depart- 
ment of Energy. The Southeastern Universities Research 
Association (SURA) operates the Thomas Jefferson Na- 
tional Accelerator Facility for the United States Depart- 
ment of Energy under contract DE-AC05-84150. We ac- 
knowledge additional research grants from the U.S. Na- 
tional Science Foundation, the Natural Sciences and En- 
gineering Research Council of Canada (NSERC), NATO, 
and FOM (Netherlands). 



[1 
[2 

[3 

[4 

[5 
[6 

[7] 
[8] 
[9] 

[10 
[11 

[12 
[13 

[14 

[15 
[16 
[17 
[18 



G. R. Farrar and D.R. Jackson, Phys. Rev. Lett. 43, 246 
(1979). 

N. Isgur and C.H. Llewellyn-Smith, Phys. Rev. Lett. 52, 
1080 (1984); N. Isgur and C.H. Llewellyn-Smith, Nucl. 
Phys. B317, 526, (1989). 

S.R. Amendolia et al., Phys. Lett. B138, 454 (1984); 

Nucl. Phys. B277, 168 (1986). 

C.J. Bebek et al., Phys. Rev. D13, 25 (1976). 

C.J. Bebek et al., Phys. Rev. D17, 1693 (1978). 

P. Brauel et al., Phys. Lett. B69, 253 (1977); P. Brauel 

et al., Z. Phys. C3, 101 (1979). 

H. Ackermann et al., Nucl. Phys. B137, 294 (1978). 
J. Volmer et al., Phys. Rev. Lett 86, 1713 (2001). 

V. Tadevosyan et al., submitted to Phys. Rev. C as Rapid 
Communication, (2006). 

T. Horn, Ph.D. thesis, University of Maryland (2006). 
R. Asaturyan et al., Nucl. Instrum. Meth.A548, 364-374 
(2005). 

J. Arrington, Phys. Rev. C69, 022201 (2004). 
M. Vanderhaeghen, M. Guidal and J.-M. Laget, Phys. 
Rev. C57, 1454 (1998); Nucl. Phys. A627 645 (1997). 
V.A. Nesterenko and A.V. Radyushkin, Phys. Lett. 
B115, 410 (1982). 

B.V. Geshkenbein, Phys. Rev. D61, 033009 (2000). 
P. Maris and CD. Roberts, Phys. Rev. C58, 3659 (1998). 
P. Maris and P.C. Tandy, Phys. Rev. C62, 204 (2000). 
A. P. Bakulev et al., Phys. Rev. D70, 033014 (2004). 



